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Abstract 

The dynamic process for the two dimensional three state Potts model 
in the critical domain is simulated by the Monte Carlo method. It is shown 
that the critical point can rigorously be located from the universal short- 
time behaviour. This makes it possible to investigate critical dynamics 
independently of the equilibrium state. From the power law behaviour of 
the magnetization the exponents f3/{vz) and l/(yz) are determined. 



1 Introduction 



Recently more and more attention is drawn to the universal behaviour of the 
critical dynamic process in the macroscopic short-time regime, for the critical 
dynamic process of model A for the 0(N) vector model starting from a random 
state with a small but non-zero initial magnetization it has been argued with 
e-expansion up to two loop that the magnetization undergoes a critical initial 
increase §. Such a scenario has been simulated for the Ising model and the 
Potts model by the Monte Carlo method |^|, [| , and the new dynamic exponent 
9 which governs the initial universal behaviour has directly been determined. 
The exponent 9 may also be estimated indirectly from the power law decay 
of the auto-correlation j|, [||, and with the damage spreading technique 
Universality |?], |j| as well as a generalized finite size scaling fa, |], |l0) has also been 
confirmed. Such an investigation of the short-time critical dynamics not only 
enlarges our fundamental knowledge on critical phenomena but also provides 
possible improvement on the numerical determination of the dynamic and the 
static exponents. The generalization and application to complex systems are 
also undergoing § 0,0,0. 

Of special interest is here a previous paper by the authors where the short- 
time dynamics for the Potts model at the critical point was simulated with the 
Monte Carlo Heat-bath algorithm ||. We have concentrated ourselves on the 
critical initial increase of the magnetization as well as the power law behaviour 
of the auto-correlation and the second moment of the magnetization. The new 
dynamic exponent 9 has directly been measured and the traditional dynamic 
exponent z as well as f3/v were also obtained. The results relate the direct 
and indirect measurement of 9 to each other and provide confident support for 
the dynamic scaling in the short-time regime. Especially the measured dynamic 
exponent z is quite rigorous. Compared with the distributed values of z between 
z = 2.2 and z — 2.7 from different simulations [jl4[, [l5|, [if], [l7|, the exponent z 



we obtained supports the small one 14, [L9 



However some important problems remain open in the numerical simulation 
of the short-time dynamics. First, the critical point was always input in the 
numerical investigation of the short-time dynamics. The critical point is usually 
defined and measured in the equilibrium state of the systems.. It is be very 
interesting to know whether one can define and locate the critical point from 
the short-time dynamics itself or not. This would allow a complete investigation 
of the short-time dynamics independent of the equilibrium state. On the other 
hand, for the critical system (with infinite volume) exactly at the critical point 
the dynamic process will, in principle, never reach the equilibrium state due to 
the critical slowing down, at least for the local dynamics. Therefore it is also 
important from the practical point of view that one can determine the critical 
point from the dynamics itself. Furthermore one may imagine that the critical 
point for the dynamic system could be different from that for the equilibrium 
state. This may conceptually broaden the scope of the critical phenomena. 
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Actually some discussion of this kind has been made with respect to the damage 
spreading |fl8|| . One purpose of this letter is to clarify how to locate rigorously 
the critical point from the short-time dynamics. 

Secondly, how to determine the static exponent (3 and v from the short- 
time dynamics in a rigorous way still requires further investigation. It has been 
pointed out that the determination of the static exponents from the short-time 
dynamics is quite promising since one does not enter the long-time regime of the 
dynamic evolution where the critical slowing down is severe I], . It has been 
demonstrated for the Ising model that the static exponents (3 and v may be 
estimated from finite size scaling. One can easily realize that the measurement 
of the static exponents from the dynamics depends essentially on the quality 
of the measurement of the dynamic exponent z. In some cases a small relative 
error for z may induce big errors for the static exponents. To determine the 
dynamic exponent z independently with the time-dependent Binder cumulant in 
principle works well. But it happens sometimes that the cumulant is somewhat 
fluctuating since it is constructed from higher moments of the magnetization. 

Therefore another purpose of this Letter is to provide an alternative way 
to determine the static exponents (3 and v from the power law behaviour of 
the observables in the beginning of the time evolution of the dynamic process. 
Compared with the method based on finite size scaling, the measurement of the 
exponents from the power law behaviour of the observables is more direct. The 
exponents may be extracted from a single lattice and the finite size effect may 
easily be under control. Actually in reference H the static exponent f3/v has 
been obtained from the power law increase of the second moment. However the 
result is still rough. What one practically measures from the second moment 
of the magnetization is id — 2(3/v)/z. As it was pointed out the value of 
(3 jv is much smaller than d and z, and therefore a small relative error for z or 
(d — 2j3jv)jz will induce a big error for 2(3 jv. We should search for a better 
way. 

Let us start with the determination of the critical point and the exponent 
\j(yz). Following the consideration in we first study the critical dynamic 
process of model A for the two dimensional three state Potts model starting 
from a random initial state with small but non-zero initial magnetization. The 
Hamiltonian is given by 

H = JJ2 <W °i = -1.0,1. (1) 

<i,j> 

According to the argument of Janssen, Schaub and Schmittmann obtained with 
the renormalization group method, one may expect a generalized scaling relation 

M^(t,T,m ) = b- k ' 3 / u M^(b- z t,b 1 / u T,b Xo m ). (2) 

to be valid also for the Potts model, for the magnetization (k = 1), taking 
b = t x l z and assuming t x °/ z mo is small enough, we can expand with respect to 
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t x °/ z mo and obtain 

M(t, t, mo) = m t e F{t 1 /^T) (3) 

where 9 = (xq — (3/v)/z. It is clear that for r = we observe the so-called 
critical initial (power law) increase of the magnetization, for r ^ 0, but still 
in the critical domain the time evolution of the magnetization will be modified 
by the scaling function Fii^/^^r). Such a fact can be used for the determi- 
nation of the critical point. In Fig. 1 the time evolution of the magnetization 
for three different coupling J = 1.0025, J = 1.0050, and J = 1.0075, with ini- 
tial magnetization mo = 0.04 and lattice size L = 72 are displayed with solid 
lines in double-log scale. The coupling for the line in the centre, J = 1.0050, 
corresponds to the theoretical critical point for the infinite system. The total 
number of samples with independent initial configurations for each coupling J 
is 600 000. By a linear extrapolation in J between the simulated curves one 
can find a particular coupling which gives the best power law behaviour for the 
magnetization. This is the "critical point" for the finite system. In principle 
there can also be some effect from the finite initial magnetization mo- The dots 
in Fig. 1 represent the time evolution of the magnetization at such a critical 
point J c = 1.0055(8). It is remarkable that the critical point located in this 
way is so close to the exact one. To estimate the critical point, we have lin- 
early interpolated the magnetization M(t, J) between the measured curves for 
J = 1.0025 and J = 1.0075. At the time t, the slope is measured within a time 
interval [t,t + 15]. Theoretically, at the critical point J c , the slope 9(t,J c ) is 
constant with respect to t. Practically the critical point is determined by search- 
ing for the minimum square deviation of the slope 9(t, J) from a constant. In 
Fig. 2 this deviation is shown when the coupling is varying around the critical 
value. The clear minimum is a confident indication for the critical point. It has 
been pointed out that in the measurement of the exponent 9, the microscopic 
time scale t m ic for the Heat-Bath algorithm is ignorably small. Therefore in 
the above estimate of the critical point, the time interval is taken to be [1, 100]. 
Actually the time interval may be taken even smaller, e.g. [1, 50]. This is really 
the physics from the short-time regime of the dynamic evolution. Furthermore, 
examining the result for smaller lattice size L — 32 and somewhat bigger initial 
magnetization mo = 0.08, one can realize that the effects of finite size and mo 
for L = 72 and mo = 0.04 are already very small, at least in the determination 
of the critical point. This is an advantage for determining the critical point 
from the short-time dynamics. 

Theoretically the critical point could be located from short-time behaviour 
of any observable in the critical domain. However practically many of them are 
not very sensitive to the coupling. Also to determine the critical point efficiently 
the dependence of the chosen observable on the critical exponents should be as 
small as possible. For example, the second moment of the magnetization and 
the auto-correlation are not good choices for the determination of the critical 
point. In reference it was mentioned that one might be able to locate the 
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critical point from the time-dependent cumulant. But this also appears to have 
some complication since higher moments are involved. 

To determine the exponent l/{yz), we differentiate In M(t, r, mo) with re- 
spect to r at the critical point, and for small enough t x °/ z m,Q 

d T In M(t,T,m )\ T=0 = t 1 '^ d T In F(t')| t <=o (4) 

Therefore d T In M(t, r, TOo)|t=q also shows a power law increase. In the nu- 
merical simulation we have approximated the differentiation by the difference 
between J = 1.0025 and J = 1.0075. This derivative is plotted in Fig. 3 in 
double-log scale. From the slope one can estimate the critical exponent l/(vz). 
Using the exponent measured from the auto-correlation as input, one obtains 
the static exponent 1/v. The result together with the exponent z and the mea- 
sured critical point is given in table 1. Due to some possible effect of t m i c , we 
have made the measurement within the time interval [20, 100]. Compared with 
the exact value 1/v = 1.2 the measured value is slightly bigger. This comes 
probably from the approximation of the differentiation by the difference. To 
improve the quality of 1/v we might represent alternatively the derivative by 
higher moments Jl9fl . Here we should point out that the slight difference be- 
tween the values for z in table 1 and that in reference || comes from some more 
serious consideration of the effect of t m i C and the finite lattice size [|| . 



Jc 


z 


(3/v 


1/v 


1.0055(8) 


2.196(8) 


0.1330(3) 


1.24(3) 



Table I: Results for J c , /3/u and 1/v obtained in this paper. The value for z used as an 
input has been taken from the previous investigation. 

Now we turn to the determination of the static exponent (5/(yz). It has 
been discussed before that it is not the best way to measure it from the second 
moment. Due to the very short spatial correlation length in the very beginning 
of the time evolution of the dynamic process the second moment behaves as 
M^ 2 \t,L) ~ L~ d . Therefore obtaining a good power law behaviour for the 
second moment is not so easy, especially for relative big lattices in the short-time 
regime. Fortunately recent investigations of the dynamic Ising model |2(], |Io| 
shows that from the power law decay of the magnetization 

M(t) ~ i~« (5) 

in the dynamic process starting from a completely ordered state, i.e. mo = f , it 
is possible to estimate the exponent (3/ {vz). With z in hand, we can obtain the 
static exponent (3/ v. Here we have simulated such a dynamic process for the 
Potts model. We have pointed out before that the critical point we estimated is 
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very close to the exact one. Therefore in the simulation here we take the exact 
value for the critical point for simplicity. The results support that universality 
and scaling appear to be valid also in a quite early stage of the time evolution 
and one can really obtain the exponent [3/{vz) quite accurately in the short- 
time dynamics. In Fig. 4 the time evolution of the magnetization is displayed 
for lattice sizes L — 72 and L = 144 in double-log scale. The total number of 
samples for independent initial configurations is 200000 for L — 72 and 80000 
for L — 144. One can hardly see any difference between both curves, i.e., there 
is already an invisible finite size effect for lattice size L — 72. From the slope 
of the curves one can estimate the exponent f3/(vz). Detailed analysis reveals 
that in this case t mic ~ 50, which is bigger than that in the measurement of 
9. However the behaviour becomes very stable after t mic . Measuring the slope 
in a time interval [100, 250] and using the value of z as input, one obtains the 
static exponent (3/v = 0.1326(6) for L = 72 and (3/v = 0.1333(2) for L = 144. 
The errors are estimated by dividing the samples into two or six groups. The 
average value is given in table 1. The static exponent f3 jv obtained in this way 
is very rigorous. In general the quality of the measurement for the first moment 
is better than that of the second moment. Besides the magnetization in this 
process is strongly self-averaging and therefore one can easily simulate bigger 
lattices. Furthermore the error of the input z will now have almost no effect on 
the quality of the exponent (3/v since the value of (3/(vz) is much smaller than 
that of z. This situation is very different from that for the measurement of the 
second moment. 

In conclusion, we have numerically simulated the critical dynamic process of 
model A for the two dimensional three state Potts model starting from both a 
random state with small initial magnetization and a completely ordered state. 
From the universal power law behaviour of the magnetization in the short-time 
regime of the dynamic evolution, the critical point and the static exponents 
are obtained. The results provide further confirmation for the dynamic scaling 
in the short-time dynamics, including the extension to the dynamic process 
starting from a completely ordered state. Especially the measurement of the 
critical point and the exponent (3/v is very rigorous. This makes it possible to 
study the short-time dynamics completely independent of the equilibrium state. 
The determination of the static properties from the short-time dynamics is not 
only conceptually interesting but also practically applicable. The application of 
this technique to the unsolved critical systems is being carried out. 
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Figure 2: Deviation from the power law behaviour when J is varied around the critical 
point. 
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Figure 3: The derivative d T In M(t, t, mo)| r =0 plotted versus time in a double-log scale. 
In order to avoid possible effects of t m j c , the slope has been calculated from the time interval 
t = [20, 100]. The dotted curve represents a straight line with the measured slope. 
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Figure 4: The magnetization M(t) versus time in a double-log scale. The solid line and 
the crosses on that line represent the results for L = 72 and L = 144. The slope has been 
obtained from the time interval t = [100, 250] in order to avoid effects from t m i c . The dotted 
curve is the straight line with the measured slope. 
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